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LOWER  3GUNDS  ON  3 EARING  ACCURACY 
FOR  CLOSELY  SPACED  SOURCES 


Abstract 


Resolution  of  multiple  sources  that  are  closely  spaced  in 
bearing  is  a  topic  of  considerable  interest  in  passive  localization  app¬ 
lications.  Many  of  the  available  studies  have  been  devoted  to  the  im¬ 
plementation  cf  algorithms  with  less  emphasis  on  the  absolute  lower 
bounds  on  attainable  accuracy.  ''This  report  presents  Cramer-Rao  lower 
bounds  on  bearing  estimation  error  for  two  closely  spaced  incoherent 
sources  which  radiate  spectrally  overlapping  broadband  Gaussian  noise. 
The  radiated  signals  are  assumed  to  occupy  the  same  frequency  band  which 
is  taken  be  to  be  lcwpass  and  rectangular  in  shape.  ^  The  signals  have 
(possibly  different)  unknown  spectrum  levels.  In  the  most  general  prob- 
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source  bearing  is  unavailable,  estimation  accuracy  is  shown  to  degrade  in 
proportion  to  the  second  power  of  the  inverse  angular  spacing,  resulting 
in  large  degradations  when  compared  to  the  accuracy  attainable  in  a  sin¬ 
gle  source  setting.  There  is  evidence  that  in  certain  situations,  a  form 
of  apricri  knowledge  concerning  signal  power  levels  can  be  used  to  signi¬ 
ficant  advantage  in  joint  bearing  estimation. 

The  detectability  of  a  second  source  from  a  signal  field  con? 
taining  two  closely  spaced  sources  that  -  overlapping/ in  both  time  and 
frequency  depends  on  the  resolving  power  of  the  array  and  is  expected  to 
be  poor  at  low  signal- to-noise  ratios  and  short  observation  times  when 
the  angular  spacing  is  small.  — .  The  problem  is  addressed  by  considering 
the  following.'  binary  hypothesis  test  : 

Hq :  There  are  two  closely  spaced  sources  present,  which  radiate 
power  SNRi  and  SNR2 ,  respectively; 

Hi:  There  is  one  source  present  radiating  the  total  received  po¬ 
wer  SNRi  +  SNR  2 • 

where  upper  bounds  on  the  false  alarm  and  miss  probabilities  are  comput¬ 
ed.  The  results  indicate  that  below  a  well  defined  threshold,  defined  by 
a  combination  of  the  two  signal- to-noise  ratios,  observation  time  and 


angular  spacing,  the  probability  of  error  in  guessing  the  wrong  hypothes¬ 
is  tends  to  one  half.  The  reachability  of  the  earlier  computed  Cramer- 
Rao  lower  bounds  for  this  regime  is  questionable.  ^ 
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Introduction 


Joint  bearing  estimation  cf  multiple  radiating  sources  is  a 
problem  which  has  generated  considerable  interest  in  the  signal 
processing  community  in  recent  years.  A  large  majority  cf  the  studies 
have  been  concerned  with  the  implementation  of  so-called  high  resolution 
algorithms,  which  have  been  known  to  be  successful  in  resolving  multiple 
closely  spaced  sources  [l]-[5j.  There  have  been  fewer  general  studies 
concerned  with  absolute  lower  bounds  on  accuracy.  Without  relevant 
performance  figures  to  serve  as  benchmarks  against  which  the  accuracy  of 
practical  algorithms  may  be  compared,  the  degree  of  optimality  for 
different  approaches  is  difficult  to  assess.  In  addition,  nest  of  the 
studies  (both  analytical  and  algorithm-based)  have  treated  the  narrowband 
problem  exclusively  [3]  -[5]. 
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er  bounds  on  bearing  error  in  a  multiple  source  setti 
ed  by  other  authors  and  most  notable  of  these  are  re 
Reference  [5]  considers  Cramer-Rao  lower  bounds  on 
11  as  the  structure  of  maximum  likelihood  bearing  es 
ases  of  both  closely  spaced  and  well  separated  source 
difference  is  smaller  than  a  beanwidth  over  the 
frequency  band,  bearing  accuracy  associated  with 
shown  to  degrade  seriously  as  result  of  the  presence 
signal  at  unknown  bearing.  The  author  demonstrates  t 
s  for  both  broadband  and  narrowband  signals.  The 
derabie  insight  into  the  structure  of  the  likelihood 
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spectral  properties  to  determine  bearing,  and 
spectral  properties  are  unknown  apriori,  a  reliabl 
son  of  the  lower  bounds  with  practical  estimation  per 
to  require  extension  of  the  unknown  parameter  set  i 
:ations  to  include  spectral  parameters. 
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Reference  [7]  presents  analytical  forms  for  the  complete  Fisher 
matrix  corresponding  to  estimation  of  both  bearing  and  power  levels  for 
two  incoherent  narrowband,  sources.  The  authors  point  out  that  for 
bearing  differences  which  exceed  a  beamwidth,  the  impact  of  uncertainty 
in  power  levels  on  bearing  accuracy  could  be  minimized  with  the  use  of 
arrays  which  exhibit  smooth  beampattern  values  in  the  neighborhood  of  an 
interference.  The  discussion  was  confined  to  source  separations  that 
exceed  a  beamwidth,  leaving  open  the  question  of  estimation  accuracy 
degradations  for  sources  in  a  closely  spaced  configuration. 


Ng  in  reference  [8]  computes  analytical  forms  of  the  Fisher  In¬ 
formation  matrix  for  a  two  receiver  array,  leaving  open  the  question  of 
performance  for  larger  arrays,  as  well  as  the  performance  degradations 
resulting  from  unknown  signal  levels.  Boehme  in  reference  [9] 
demonstrates  with  extensive  simulation  results  the  statistical  coupling 
of  bearing,  range  and  signal  level  estimates  in  a  multiple  source 
setting.  His  study  considered  both  well  separated  and  closely  spaced 
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ces  but  presented  no  performance  figures. 

This  study  considers  the  bearing  estimation  accuracy  of  tw< 
ely  spaced  sources  under  conditions  in  which  prior  knowledge  of  th< 
ated  signal  powers  is  not  necessarily  available.  The  source  signals 
modelled  as  spectrally  and  temporally  overlapping.  The  stud] 
s tigat.es  the  performance  of  line  arrays  of  equally  spaced  sensors 
r  ccnditicns  of  practical  interest.  The  power  ratio  cf  the  tw< 
ce  is  demonstrated  to  play  a  major  role  in  determining  accuracy.  Th) 
ussion  is  largely  focussed  on  the  situation  in  which  the  angulaj 
erence  between  the  sources  is  extremely  small. 

At  low  signal  to  noise  ratios  and  short  observation  times,  th) 
ity  cf  the  array  to  distinguish  the  two  sources  from  a  single  sourci 
expected  to  deteriorate,  and  this  suggests  that  without  apriori 
'ledge  of  the  number  of  sources  present,  the  Cramer-Rao  lower  bounds 
simply  unreachable.  The  final  portion  of  the  study  addresses  thl 
liability  of  the  two  sources  when  the  angular  difference  becomej 


small. 

The  report  is  organized  as  follows.  In  section  1,  the  genera! 
theory  is  presented  which  includes  analytical  forms  of  the  lower  bound: 
on  source  bearings.  Section  2  discusses  features  of  the  bounds  when  th 
signal  levels  are  known  apricri.  The  relative  signal  to  noise  ratios  o: 
the  individual  sources  is  shewn  to  play  an  important  role  in  determinim 
overall  accuracy.  Unknown  signal  power  levels  are  introduced  in  sectio: 
3  where  accuracy  reductions  in  the  performance  computed  in  section  2  ar 
shown  to  be  quite  large.  Section  4  considers  the  resolvability  of  th 
sources.  A  summary  of  the  important  results  is  given  in  section  5. 
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incoherent  radiating  sources  (Figure  1)  . 


source 


The  signal  received  by  any  given  sensor  in  the  array  i 
of  the  two  source  signals  in  addition  to  a  noise  signa 


x$.(t)  =  s1(t-ti)  +  S2 ( t-Zi)  +  ni(t) 


i  =  x, 


It  is  assumed  that  the  signals  s^ft)  and  S2(t)  are  mutually  incoheren 
but  individually  (perfectly)  coherent  over  the  receiving  array,  while  th 
noise  signals  are  statistically  independent  from  sensor  to  sensor.  '  Th 
variable  tj_  denotes  the  propagation  delay  for  source  signal  1  to  trave 
to  sensor  i.  Similarly  z^  refers  to  the  delay  associated  with  sourc 
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Figure  1 


The  sensor  spacing,  d,  is  uniform  over  the  array.  Since  there  are  a 
total  of  M  receivers  in  the  array,  the  array  has  length  (ti-l)d  i  L.  The 
angular  difference  between  the  sources  measured  with  respect  to  the  array 
axis  at  the  midpoint  is  La. 

Based  on  signal  observations  from  the  sensors,  we  wish  to  fin 
the  lower  bound  on  estimates  of  and  «2*  If  the  spectral  properties  o 
the  signal  and  noise  processes  are  known  apriori,  the  parameters  o 
interest  include  only  bearings  and  the  relevant  parameter  vector  9.  is 
given  by  : 


£  =  [  ai,  02]  T 


The  inverse  of  the  well  known  Fisher  Information  matrix  [11]  yields  the 
error  variances  of  estimates  for  parameters  in  6.  : 
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Cov  {  5if5j  )  >  J*1  ji.  (3) 

If  the  spectral  parameters  are  unknown,  e_  will  be  suitably  enlarged.  In 
the  first  part  of  this  discussion,  the  effects  of  uncertain  spectral 
parameters  will  be  ignored. 

The  following  assumptions  are  used  in  the  analysis  : 

1.  The  signals  S]_(t)  and  S2(t)  are  sample  functions  from 
statistically  independent  stationary  Gaussian  processes  with  zero  means. 

processes  have  lowpass,  rectangular  spectra  of  equal  width,  W ,  but 
with  possibly  different  heights. 

2.  The  noise  signals  appearing  at  the  various  receivers 
are  sample  functions  from  statistically  independent,  stationary 
Gaussian  processes  with  zero  means  and  identical  soectral  orooerties 
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The  observed  data  is  represented  by  the  vector  of  Fourier  ccef- 
iiciar.ts  of  the  signals  received  by  the  various  sensors  concatenated: 


where 


X  =  [  X(o)]>)^',  X(w2^' 


X(ujc)  =  [  Xi  (ufc)  ,X2  («fc)  » 


,  X  (  ujj)  ^  ]  T 


,  Xm(«Jc)  3T 


Xi («k)  = 


Xi(t)  e 


The  probability  density  function  of  the  data  vector,  X,  condit¬ 
ioned  on  8.,  is  complex  and  Gaussian  with  zero  mean.  If  the  time- 
bandwidth  product  TW  is  much  larger  than  one,  the  data  Fourier  co¬ 
efficients  associated  with  different  frequencies  are  approximately  un¬ 
correlated.  Under  these  cond-itions ,  the  pdf  of  the  data  vector  X  has 
the  form: 
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p(X/9)  = 


C  dec  iiXk  ]_1  exp  (  -X;<  K*"1^ 
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(4) 


the  covariance  matrix  of  the  sensor  data  at  frequency  4^ .  Define 
gnal  steering  vectors  yi_  and  22  as  follows  : 

VI  -  [  . ]T  (5) 

5t2  *  C  e-j“^zl . ]T  '  (5) 

n  has  the  form  : 


Xjj  =  S  y  (  )  Vi_yi_ 


S2  (  «}<)  X.2X.2  +  N(  a^)! 


(7) 


I  is  the  identity  matrix.  S^U^),  S2(«]c)  and  NU^)  are  the  power 
a  of  S]_(t),  s2(t)  and  n^(t),  respectively. 


The  elements  of  J  can  be  obtained 
ig  well  known  results  (11)  : 


from  equation  (2) ,  using  the 
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Kfc  can  be  inverted  with  the  use  of  Bartlett's  Identity: 
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The  k  -  subscripts  have  been  omitted  in  the  quantities  Pj_  and  P2  in  (10) 
to  simplify  subsequent  algebraic  expressions  : 


p_  _  Si  ( )  /M  (  a»y  ) 
^  1  +  MSi_(«]^) 


p2  = 


S2  (av)  /N  ( ) 
1  +  MS 2 ( ) 


izi"z2  i 2  =  c< 

i  /  j  =  1 
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Ml;  (  ti*  t  i  ~  Zi  *  Z  i  ) 


(12a) 


(12b) 


Note  the  dependence  of  the  function  } v-|  ~vo  j  on  the  angular  frequency  a 
and  the  differential  bearing  and  range  between  sources  which  is  expressed 
by  the  quantity 

(  tj_  ~  t  j  —  2  ~  Zj  )  . 

Calculations  carried  out  in  Appendix  1  yield  the  following  exp¬ 
ressions  for  elements  of  the  Fisher  Information  matrix  : 
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J11  =  ^  Gk  •  P1  Si  Uk)  /N(uk)  {  31  "  p2i’l  j  +  pl2p2 
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k  =  1  »*  J 

The  variables  >?i ,  n2 ' 21  *  z2 » *1  *  s2  are  defined  as  follows  : 
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3j_  is  recognized  as  the  familiar  pattern,  function  derivative. 
From  an  analytical  point  of  view,  is  the  only  geometrical  quantity 
that  is  important  in  determining  the  accuracy  of  a  location  parameter  for 
a  single  source.  One  can  show  this  by  setting  S2(«^)  equal  to  zero  over 
the  frequency  band  in  (13)  .  It  then  follows’  that  ■*  (  1  + 
MSi  (  /N  ( 4)]^)  )  and  one  obtains  for  : 


Var  (ai  )  =: 


(Si  (ait)  /NUy. )  )  2 

+  MSi_  (  u^)  /N  (  ufc ) 


(19a) 


It  can  be  shown  (  see  Appendix  2)  that  when  the  sensors  are  arranged  in  a 
line  array  configuration  of  total  length  L,  (19a)  is  equivalent  to: 
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Var  ( Si  ) 


%  f  \  (Si  (ut>)  /N(«v)  ) 

~  \  1  +  MSI.  (aj^) /Nl«> 


k  =  1 


r  '  >  }  ’l 


(19b) 


where  Xjj  is  the  signal  wavelength  at  frequency  k.  The  above  expression 
for  the  variance  of  bearing  error  indicates  that  the  important  geomet¬ 
rical  quantity  determining  accuracy  is  the  ratio  Lsin(a^) /Xj*.  LsinCa^) 
is  the  array  baseline  (  length  of  the  array  component  in  the  direction 
orthogonal  to  the  line  of  sight  to  the  source  ) .  Accurate  bearing  esti¬ 
mates  will  therefore  result  when  the  baseline  exceeds  a  signal  period. 
Since  the  Cramer-Rao  lower  bound  is  a  local  error  bound,  it  will  not 
yield  information  about  estimation  accuracies  resulting  from  ambiguous 
estimates  (  noise  generated  peaks  in  the  likelihood  function  which  *  occur 
far  away  from  the  true  parameter  value,  or  multiple  solutions  to  the 
likelihood  equation  resulting  from  ambiguous  delay  estimates  for’  narrow- 
band  signals  ) .  If  the  possibility  of  ambiguity  error  can  be  safely  dis¬ 
counted,  then  equation  (19b)  states  that  reductions  in  the  Lsinfai/Xfc) 
ratio  must  be  offset  by  simultaneous  increases  in  the  signal  to  noise 
ratio  or  observation  time,  if  operation  within  reasonable  accuracy  limits 
is  to  be  acheived.  The  remainder  of  this  report  will  be  confined  to  a 
discussion  of  local  errors  only. 

/ 

Returning  now  to  equations  (13)  -(18),  it  is  clear  that  the  remaining 
quantities  appearing  in  the  Fisher  matrix  describe  the  interactions  of 
the  two  sources  resulting  from  effects  of  interference.  Of  prominence 
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The  be amp at tarn  is  defined  in  equation  (12b)  and  will  be  shown  to  depend 
on  the  signal  frequency,  the  array  configuration  and  the  angles  of 
orientation  of  the  two  sources.  The  array  bear.pattarn  is  the  measured 
response  at  a  particular  frequency  for  an  array  which  has  been  focussed 
urce  1  but  which  receives  signal  energy  from  the  direction  of  an 
erer.ce  (source  2).  At  frequencies  for  which  the  bear.pa t tern  is 
the  sources  appear  to  radiace  fror.  nearly  the  sar.e  bearing  angle, 
equencies  for  which  the  bear.pattarn  is  low  the  sources  are  easily 
guishable.  Of  particular  interest  to  us  is  the  case  in  which  the 
s  appear  to  er.ar.ate  fror.  nearby  bearings  for  all  signal  frequencies 
signal  band. 

It  can  be  easily  shown  with  a  few  seeps  of  algebra  that  if  the 
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le,  che  function  C-^  varies 
with  frequency  and  the  resulting  k-sums  r.ust  be  evaluated  with  care.  It 
will  be  convenient  to  consider  bearing  accuracy  under  specific  conditions 
of  interest  and  for  these  cases  the  behavior  of  the  k-sur.s  can  be  rade 
sir.ple.  The  next  section  considers  the  probler.  of  closely  spaced 
sources . 

2.  Bearing  Accuracy  for  Closely  Spaced  Sources 

Consider  now  the  problem  of  closely  spaced  sources,  where  oj_ 

*2  is  sufficiently  small  so  that 

2  2 

|vifv2!  s  M  .  (20) 

If  the  source  signals  impinge  on  the  array  in  the  form  of  plane  waves, 
elements  of  the  Fisher  matrix  can  be  simplified  considerably.  Choosing 
the  array  midpoint  as  the  coordinate  system  origin,  for  plane  wave 
arrivals  the  travel  times  are  given  by  : 

tj_  =  idcosa^/c  ;  =  idcosa2/c  ;  i  =  -  (M-l) /2 ,  .  .  .  ,  (M-l) /2  (21a) 

ti  -  Zi  =  id(cosai  -  cosatj)/c  (21b) 

The  beampattern  becomes  : 


|Vi*V2  ;  =  y  cos  j  (ujv/c)  (l  -  j)d(cosay  ~  ccsay)  f 

^  j 

.  .  M-l 
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Suppose  jay  -  122!  <<  2s  : 


ty  -  zy  S  idAasir.a ;  Act  =  ay-a2  i 


a  =  ( ay  +  a 2 ) / 2 ; 


(23a) 


Under  these  conditions  the  cosine  argument  of  the  function  in  (22)  wil 
vary  negligibly  over  the  i  -  index  provided  AaLsin(a)W/c  <<  2s.  In  sucu 


a  case 


2  ■  2 

1-  Ivitvol/M  <<  1 


(23b) 


Note  that  the  condition  (23b)  would  be  satisfied  in  general  when 


(ufcd/c)  (i-j)  (cosay-COSa2)  1  <<  2s  v  k 


In  addition  to  the  case  ay  s  a2 .  the  criterion  of  equation  (24)  can  also 
be  sec  whenever  L  <<  \~yn.  We  assure  that  this  is  not  the  case,  that  is, 
we  wish  to  be  in  an  estimation  regime  in  which  accurate  estimates  of  ay 
or  a2  are  possible  when  the  measurements  are  obtained  without  the  pres¬ 
ence  of  the  second  signal  arrival.  When  L  <<  \ain  t*le  array  simply  has 
no  means  of  determining  bearing  even  for  a  single  source  since  differen¬ 
tial  delays  which  are  fractions  of  a  signal  period  cannot  be  exploited  in 
bearing  measurement. 

Consider  then  the  lower  bound  for  closely  spaced  sources,  i.e.,  equat¬ 
ion  (24)  holds  but  L  >  *nxn*  Computations  carried  out  in  Appendix  2 
yield  the  following  approximate  expressions  for  the  Fisher  matrix  ele¬ 
ments.  When  contributions  to  the  lower  bound  of  order 

2  2 
(0W  ) 

and  smaller  are  neglected  one  obtains  the  following  : 
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f  the  TW  product  is  large  compared  to  one,  the  k  -  so 
e  above  expressions  can  be  approximated  by  integrals. 
Define  the  following 
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(31) 


(32) 


In  terns  of  the  above  integrals,  the  F.I.M.  (Fisher  Information 
elements  can  be  expressed  as  follows  : 


matrix) 


Jll=  X12(MS1/N)2j  li  +  p (MS2/N) 2  I 3  )  (33) 

J 

J12  -  XiX2M2S1S2/N2(  It  -  (6/5)  0I2  -  /5M2S!S2/N21 3  j  (34) 

J22  *  X22(MS2/N)2i  II  +  P (MS^/N) 2  I 3  )  (35) 

j 

Consider  equations  (33)  —  (35)  ,  which  are  explicit  functions  of  I-t,  I7  and 
I3.  If  * 

2  2 

1  +  MS/N  >>  M2S!S2/N2  ^  S  M2S1S2/N  (  1  -  5|_  +  .  .  .  )  <36) 


over  the  frequency  band  then  the  denominators  of  the  integrands  are  dom¬ 
inated  by  the  constant  term,  1  +  MS/N.  Under  this  condition,  the  inte¬ 
grands  of  I]_,  I2  and  I3  will  grow  (  approximately  )  with  «2,  u4  and- •  <b°, 
respectively.  If,  on  the  other  hand 
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1  +  MS/N  <<  M2S!S2/N25  =  M2SxS2/N2  (  1  -  £|_  +  .  .  .  )  (37) 

over  some  frequency  band,  then  integrands  of  I1,I2  and  I3  will  tend  to 
vary  {  approximately  )  with  ,  «>2  and  respectively.  Note  that  there 
is  no  other  possibility  for  the  behavior  of  H^. 

Obviously  (37)  will  hold  only  at  large  signal  to  noise  ratios. 

In  fact  in  order  for  (37)  to  hold 


MSi/N  >>  Ipoif2)  ~1 


i  =  1  and  2 


1 


For  example,  if  <  .1,  then  MSj_/N  >  10  for  i  =  1,2.  On 

hand,  equation  (35)  will  hold  when 


MSi/N  <  (Jmjc2)"1 


for  i  =  1  or  2 


Define  the  variable  R  : 


M2Si So/N2 


1  +  KS/N 


We  will  consider 
extreme  values  : 


two  regir.es  of  interest  occurring  at 


<<  1 


>>  1 


where  £r.ax  is  defined  as  the  value  of  5  occurring  at  the  largest  (co 
signal  frequency,  w. 

When  3  5nax  <<  1  it  is  shown  in  Appendix  3  that  eler.ents  of 
F.I.M.  can  be  approximated  by  : 
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After  a  few  steps  of  algebra. 
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3y  daf  ir.it.icn ,  MSj/N  <<  (0W  )  sc  that  the  bearing  e 
j  is  much  larger  than  for  source  i.  In  order  to  obtain 
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(19b 


the  effects  of  the  presence  of  an  interference  on  the  abi 
mine  bearing,  it  will  be  useful  to  compute  the  fractional 
(aj_)  ,  defined  as  the  ratio  of  bearing  error  for  estimates  of 
ed  under  the  assumption  that  source  j  {  the  interference  )  is 
located  at  an  unknown  bearing,  to  the  equivalent  performance 
ir.ed  with  source  j  absent.  The  lower  bound  corresponding  to 
situation  is  given  by  expression  (19b) .  With  the  definition 
)  is  equivalent  to  : 


msigr.t 
lity  to 
error , 
a  ^  cb  — 
present 
figure 
the  lat- 
for  Xj_, 


Var  ( on /source  j  absent)  i  (MS^/N)‘‘Ii 


S2  =  0 


,  r  TW  (MS--/N)2  W2  „  2  1  "1 
"  L  2  t  1  +  MS ,  /N  3  1  j 


Therefore  we  have  : 


AVar  (Si)  i 


Var ( S  ^  / source  j  present) 
Var  ( S^ /source  j  absent) 
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For  the  low  SNP.  case  :  MS/N  <<  ( ) 


2  -1  f  2.03 
AVar  (SO  =  (MSi/N  /JW  )  *  ( 

1.66 


;  ms/n  <<  i 

2  -1 

;  1  <<  MS/N  <<  (0W  ) 
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For  MSi/N  >>  (/JW  )  ;  i  =  1  or  2  : 


AVar  (Sf)  =  1 

2  -1  (47) 

AVar  ( 5  j  )  =  5  (MS  j  /N  /JW  )  j  -  i 

Equations  (43)  and  (46)  show  that  the  value  of  the  bearing  error 
associated  with  is  insensitive  to  the  signal  to  noise  ratio  S2/N, 

since  the  total  variation  of  the  bound  over  the  summed  SNR.  is  smaller 
than  2  db.  The  resulting  bearing  error,  however,  is  larger  by  the  factor 
(/JW2)~1  when  compared  to  the  case  of  source  2  absent.  Equations  (44)  and 


(47) ,  on  t he  other  hand,  shews 
curacy  in  the  estimation  of 
(MSx/N)*0W*  >>  l  provided  that 


that  one  can  obta 
ox  if  is  suf f 
Sj  is  much  small e 


undegraded  bearing  ac- 
i an t ly  large  so  that 
than  S-;  . 


details 


The  case 
leading 


*£nax  >>  1  is  considered  next.  Appendix  4  contains  the 
to  the  following  expression  for  the  Fisher  matrix  : 
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It  is  straightforward  to  show  that 
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(43) 


(49) 


Equation  (49)  decreases  mor.otonically  as  the  ratio  Sj/Sj  tends  to  zero. 
The  fractional  error  cr  penalty  on  bearing  accuracy  resulting  from  the 
presence  of  a  second  source  in  the  estimation  of  either  ox  or  02  is  given 
by  : 


AVar(Si)  «  (JW  /12)  (  1  -  ) 

S2 


(50) 


As  Sj/Sx  tends  to  zero. 


fi\*  MSj/N 


AVar  { Sj. ) 
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(51) 


A  Var  (  a  ) 


fiV  MS  j_/'t 
12 


(52) 


Since  both,  products  or.  the  right  sides  of  (51)  and  (52)  are  larger  then 
one  by  assumption,  the  fractional  increase  in  bearing  error  for  *4  >>  1 
is  quite  large. 

Figure  2  depicts  the  fractional  error  for  bearing  a}.  with  MS^/M 
<<  (/!W2)-1  {  solid  curve  )  and  MS^/N  >>  (£w2)~1  (  dashed  curve  )  as  a 
function  of  MSs/N.  The  accuracy  reductions  are  significant  at  all  values 
of  MSj/N  for  the  solid  curve,  but  present  only  at  the  high  end  for  MS-;/N 
in  the  dashed  curve. 
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Figure  2 

I  t  is  worth  pointing  out  that  if  one  of  the  sources  happened  to 
be  at  a  known  bearing,  the  lower  bound  on  estimates  of  the  unknown  source 
bearing  has  the  following  interesting  properties.  Let  source  j  be  loca¬ 
ted  at  known  bearing  and  source  i  be  at  unknown  bearing.  If  source  j  has 
a  signal  to  noise  ratio  which  is  jauch  smaller  than  (/jw2)“1  the  degrada¬ 
tion  factor  is: 


Judging  from  (53a),  the  incremental  bearing  error  risas  linearly  with  or. 
raoio  of  ir.carferancs  to  signal  powers.  However,  if  the  SNR  of  scuroa  j 


MS-i/N 


exceeds  (jJW-)--  as  well  as  MS-/N,  the  degradation  factor  is: 


1  +  Sj/S.-  2  -1 

AVar(Si)  =  - - —  =  (MS.-/N  fiW  )  >>  1  (53b) 

MSj/N  fiV2 

2 

Since  (MSj/N)  0W  >>  1,  (53b)  is  smaller  than  (53a)  but  large  degrada¬ 

tions  in  bearing  error  remain  as  result  of  the  presence  of  source 

2  -1 

2.  If  the  signal  to  noise  ratio  for  source  i  is  large,  MS-;/N  >>  ( 0W  )  , 

2  -1 

but  MSj/N  <<  (£W  )  ,  no  degradations  in  accuracy  are  observed. 

The  relationship  between  (53a)  and  (53b)  can  be  explained  as 
follows.  Suppose  Aa  =  0  --  this  corresponds  to  (53a).  Any  finite  MSj/N 
must  be  smaller  than  prJ2  since  p  is  zero.  As  Aa  increases,  a  large  value 
of  MSj/N  exists  for  which  MSj/N  >>  (0W2)“1.  However,  as  Aa  increases 
sufficiently  the  beampattern  will  tend  to  zero  and  the  function  5  will 
tend  to  one.  At  this  point  there  is  no  accuracy  reduction  in  estimating 
<*2.  because  the  sources  have  become  well  separated.  Hence  the  transition 
from  (53a)  to  (53b)  describes  the  effect  of  increasing  separation  on 
bearing  accuracy  improvement. 
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When  both  signal  to  noise  ratios  exceed  ip W  )  one  can  compute 
analogous  degradation  factors  : 


2 

1.5  MS2/N  /JW  ;  Si  >>  Sj 
2 

.75  MSj/N  jJW  ;  Sj  >>  Si 

2  -1 

With  MSi/N  >>  (fiW  )  ,  accuracy  reductions  resulting  from  the  presence  of 

a  second  source  are  large  in  this  setting  also. 

At  low  signal  to  noise  ratios  for  the  source  at  unknown  bearing, 
increases  in  the  signal  to  noise  ratio  for  the  source  at  known  bearing 
tends  to  improve  bearing  accuracy.  At  large  signal  to  noise  ratios  for 
the  source  at  unknown  bearing,  reductions  in  accuracy  are  present  only 
when  the  SNR  of  the  other  source  is  large  as  well.  In  this  setting,  the 
optimal  condition  is  =  S2,  where  the  error  increases  monotonically  as 
the  power  ratio  deviates  from  one.  Hence,  even  when  the  bearing  of  one 
source  is  known  apriori  ,  the  ability  to  determinine  the  bearing  of  a  se¬ 
cond,  closely  spaced  source  can  be  quite  poor. 

The  results  of  this  section  indicate  that  a  potentially  useful 
regime  exists  in  the  estimation  of  a  particular  source  bearing  provided 
that  the  SNR  for  that  source  is  sufficiently  large  and  the  interfering 
source  SNR  is  sufficiently  small.  When  the  conditions  ^.fail  to  hold. 
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bearing  accuracy  is  greatly  reducad  in  the  multiple  source  setting.  In 
the  next  section  it  will  be  shown  that  unknown  spectral  levels  corr.pror.ise 
bearing  accuracy  further  and  that  under  these  conditions  one  is  always 
forced  to  operate  with  large  accuracy  reductions. 


3.  Bearing  Estir.ati 
Levels 


on  Accuracy  Without  Prior  Knowledge  of  Source  Power 


When  prior  knowledge  of  the  source  power  levels  is  unavailable 
the  relevant  parar.eter  vector  9.  becor.es  : 


1  =  l  al '  o<2  -  S1 '  s2  ]T- 
Consider  the  following  partitioned  r.atrix 
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where  is  the  Fisher  r.atrix  associated  with  bearings  only 
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The  total  bearing  error  can  be  expressed  as  follows  : 


Cov(3i,Sj)  > 
Define  the  matrix  X  : 
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Appendices  6-3  contain  details  leading  to  expressions  for  the 
F.I.M.  under  various  conditions.  It  is  useful  to  consider  properties  of 
the  J3  r.atrix  alone  before  discussing  the  problem  of  joint  bearing  and 
power  estimates.  The  following  statements  are  based  on  derivations  taken 
from  the  appendices. 


S2  =  0 
lated . 
signal 


1.  When  the  signal  field  contains  only  a  single  arrival  (e.g., 

,  estimates  of  bearing  and  signal  power  are  statistically  ur.corre- 
The  lower  bound  on  the  variance  of  signal  power  estimates  for  a 
of  strength  S  is  given  by  : 


Var  ( )  > 


TW 
2  3 


2 

(M/N) 

(1  +  MS/ N) <• 


-1 


a 


1 

TW 


N2/  M2  MS/N  <<  1 
S2  MS/N  >>  1 


(31) 
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2.  When  the  signal  field  contains  arrivals  from  two  incoherent 
sources,  bearing  and  signal  power  estimates  are  in  general  statistically 

2  -1 

correlated.  For  signal  to  noise  ratios,  MSj/N,  smaller  than  ( jJW  )  : 


Var (Si) 
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(83) 


2  -1 

The  degradation  factor  is  proportional  to  (jJW  (1  +  MS/N))  >>  1. 
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3.  If  one  signal  to  noise  ratio  exceeds  (jJW  )  (  for  example, 

source  i  )  estimates  of  signal  powers  are  statistically  uncorrelated. 
The  lower  bounds  are  given  by  : 
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If  one  defines  integrals  I4  and  I5  as  follows,  the  elements 
J2  and  J3  can  be  expressed  entirely  in  terns  of  integrals  Ii  -  I5. 


rn 

-4  4  17  I  TT 


2  w 

<t)  G  uj 

MS/N  +  M*SiS2/N*  5) 


15  4  57  f  IT 


+  MS/N  +  M*SiS2/N*  5)  *• 


With  the  following  expressions  front  Appendix  3  for  nx  and  ;j2  ; 

1/2  2  2 

Ii  =  4)  (  1  -  2/5  fin  +  .  .  .  ) 

one  can  show  that 
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MSi/H  >>  (fiW  ) 


the  lower  SNR  source  has  poorer  accuracy  in  spectrum  level  estimation 
than  the  large  SNR  source.  Both  estimation  errors  exhibit  degradations 
resulting  from  the  multiple  source  setting,  in  spite  of  the  fact  that  the 
errors  are  actually  uncorrelated. 
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4.  If  both  signal  to  noise  ratios  exceed  the  value  ( fiW  )  ,  the 

estimation  errors  of  Si_  and  Sj  are  also  statistically  uncorrelated  : 


^  t*V  ; 


;  i  »  1,2 


We  consider  now  the  joint  parameter  estimation  oroblem.  Suptose 
2-1 

KSi/N  <<  ( fiW  )  ;  i  =  1,2.  Computations  carried  out  in  Appendix  6  for 

s  of  the  X  matrix  are  as  follows  ; 
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Tha  incremental  error  due  to  uncertainty  concerning  signal  power  levels 
is  defined  as  follows  : 
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At  large  signal  to  noise  ratios,  MS/M  >>  1  : 
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The  incremental  error  due  to  unknown  power  levels  is  given  by: 
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1  +  5  S  j_S  2/  3  S 


which  grows  with  MS/N. 

At  low  signal  to  noise  ratios  the  incremental  errors  due  to  un¬ 
known  amplitude  factors  is  a  constant  value  of  approximately  5 . 4  db .  In 
tha  large  signal  to  noise  ratio  case  the  errors  are  likely  to  be  even 
larger  and  this  suggests  that  apriori  knowledge  of  signal  levels  might  be 
used  to  advantage  in  improving  overall  accuracy. 

2  -1 

Next  consider  the  following.  Suppose  MSi/N  <<  ( pw  )  but  MS2/N 
2-1 

>>  (0W  )  .  Appendix  7  computes  the  following  lower  bounds.  At  low 

signal  to  noise  ratios  for  source  1,  i.e.,  MS^/N  <<  1  : 
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Var  (a2)  *  (  J7  X22  W2  /JW2  y|  MS2/N  [ 

from  which  the  incremental  bearing  errors  can  be  obtained  : 

5  Var  ( )  =  =  5  V a r  ( a 2 ) 

**1 

At  large  signal  to  noise  ratios  for  source  1  {  MS]_/N  >>  l) 
Var  (Si)  >  |  X!2  W2  |  fiV2  (MSx/N)2  |i  } 

Var  (52)  -  {  57  x22  W2  |  *»2  M2S1S2/M2  j 

The  incremental  bearing  errors  are 
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The  effects  cf  uncertain  signal  levels  produce  iarg 
degradations,  particularly  for  5 2  which  previously  exhibited 
tal  bearing  errors  resulting  fro:?,  the  presence  of  source  1. 
Finally  we  consider  the  last  case  : 
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The  incremental  bearing  error  resulting  from  lack  cf  knowledge  of  power 
levels  is 
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2 

;  Si  <<  Sj 

The  above  degradation  factor  increases  monotonically  with  the  ratio 
Sj/Si.  It  is  also  proportional  to  {Kfi which  is  by  definition,  a 
factor  exceeding  one. 

Figure  3  depicts  the  incremental  bearing  error  due  to  unknown 
signal  levels  for  3i  as  a  function  of  MSj/N  for  MSj_/N  <<  (/Jw2)~l  .  (solid 
curve)  and  MS^/N  >>  (dashed  curve).  In  the  lower  SNR  case  for 
source  i,  the  incremental  error  increases  monotonically  with  MSj/N,  un¬ 
like  the  nearly  constant  behavior  of  the  fractional  error  for  known  sig¬ 
nal  levels  of  Figure  2.  In  the  large  SNR  case  for  source  i  the  incremen¬ 
tal  error  decreases  as  MSj/N  increases. 
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wr.ere  resolvability  of  the  sources  ar.d  reachability  of  the  Cra.ter-Rao  lo 
wer  bounds  is  questionable.  The  next  section  is  concerned  with  the  errc 
associated  with  a  binary  hypothesis  test,  where  there  is  one  source  pres 
r.der  the  first  hypothesis  and  two  pr 
bility  of  error  in  choosing  the  in 
es  are  expected  to  be  distir.guisha 
ution  of  the  sources  becor.es  unlike! 
r-Rao  lower  bounds  will  be  poor. 

robability  of  Error  in  Resolving  the 

Determination,  of  the  correct  nur.se r  or  sources  present  using 

s  study  because  the  pr 
is  reasonable  to  expect 
SNR  and  lone  observation  tire,  two  sources  ca 


nowever,  one  expects  that  the  ability  to  precict  the  existence  of  the  se¬ 
cond  nearby  source  will  deteriorate  as  the  SNR  and  observation  tire  be¬ 
come  small.  This  section  considers  this  problem  from  the  viewpoint  of 
determining  the  probability  of  error  in  choosing  the  wrong  decision  in 
the  following  binary  hypothesis  test  : 


K0  : 

Xi (t)  =  (Si+S2) 
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Hi  : 

Xj_(t)  =  (Si) 

Under 

Ki  there  are  two 
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ng  angles  ai  and 
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cated  at  bearing  ag  ,  where  ag  is  the  "weighted  average"  bearing: 
sial  +  s2a2 

a o  =  -  ;  and  =  dsin(ag)i/c. 
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;  and  vj_  =  dsin(ag)i/c 


Note  that  if  Si  =  S2,  ag  =  a,  the  true  average.  If  Si  >>  S2  so  that  S  2 
Si  then  ag  s  ay.  The  signals  si(t),  S2(t)  and  s(t)  as  assumed  to  have  i- 
dentical  lowpass  spectra  of  width  W  and  with  unity  height. 

The  following  material  is  contained  in  reference  [7],  pp .  - 
where  lower  bounds  on  the  false  alarm  and  miss  probabilities  in  a  binary 
hypothesis  test  are  presented. 

Define  the  probability  density  functions  of  the  received  data 
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follow! 


Based  on  calculations  carried  out  in  Appendix  3  one  obtains  tie 
ng  expression  for  ji(s).  Define  : 
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Equation  (107)  is  a  function  of  the  two  steering  vectors  vi_  and 
V2  as  well  as  the  third  vector  vq'  w^ich  is  still  perfectly  general. 
Next  let  vq  represent  the  steering  vector  of  delays  for  a  signal  radiated 
at  bearing  aQ  : 

Si«i  +  S2«i2 


After  simple  computations  carried  out  in' Appendix  8  one  can  show  that 
ji(s)  takes  the  form  : 
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0  depends 


n  order  to  obtain  the  solution ,  it  is  useful  to  consider  the 
he  quantity  ?.  y«~ .  Suppose  the  angle  difference,  ii,  upon  which 
happened  to  be  extremely  snail.  For  any  finite  signal  to  noise 


ratios  (and  hence  f ini oe  value  of  ?)  there  exists  a  value  of  Act  su 


will  be  smaller  than  one.  It  will  become  apparent  later  that  a  threshold 
region,  above  which  the  error  probability  tends  to  zero,  will  depend  u 
or.  a  variety 'of  parameters.  A  low  probability  of  error  indicates  th 
resolution  of  the  sources  is  likely.  An  error  near  cr.e  half  indicates 
that  the  data  is  essentially  useless  in  determining  the  number  of  sources 
present.  Since  the  minimum  angular  separation  for  which  the  two  sources 
can  be  distinguished  is  of  ultimate  interest  to  this  study,  we  choose  Aa 
so  that  iy*-  is  much  smaller  than  unity.  Then  the  parameter  combinations 
necessary  to  cause  the  probability  of  error  to  vanish  can  be  determined 
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/\riirJ  =10,  a  -  t/2,  TV  =  100,  and  MSi/N  =  MS 2/N  =  .05,  resolution 
ssible  when  4a  >  .15  radians  =  8.6  degrees.  If  the  signal  to  noise 
s  increase  to  .25  each  but  with  the  rer.air.ing  parar.eters  kept  fixed, 
ir.ir.ur.  resolvable  angle  difference  is  5. 73  degrees. 

Detectability  of  two  scurc 

y  spaced  sources  that  are  overlapping  both  ter.poraily  and  spectral 
related  to  the  resolving  ability  of  the  array.  This  section  ha 
ir.ed  the  necessary  parar.eter  combinations  that  lead  to  accurate  de 
n.  Under  such  conditions,  the  two  signals  are  distinguishable  an 

r-Rao  lower  bounds  can  be  reached.  When  th 
(when  one  is  operating  belo 

wo  sources  are  essentially  indistinguishable.  Without  apriori  know 
of  the  correct  number  of  sources  present,  one  simply  cannot  hope  t 
the  Cramer-Rao  lower  bounds. 
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scur.zs  cr.  :e  =  r;r.;  accuracy  wear.  :v:  c.cse.y  spacer  ir. cc rarer. c  sources 
radiate  broadband  Gaussian  noise  and  are  received  by  an  equally  spaced 
line  array  in  the  presence  cf  spatially  wr. ice  Gaussian  noise .  The  sour¬ 
ces  are  taken  to  he  closely  spaced  whenever  chair  difference  in  bearing 
is  much  scalier  than  a  be  a. -.width  .  The  sources  are  overlapping  both  ten¬ 
or  a  1 1  y  and  spectrally.  The  spectral  shapes  of  the  signals  are  taken  to 
owpass ,  rectangular  functions  of  identical  width,  W„  but  with  possibly 
ifferent  heights.  The  uni  fern,  sensor  spacing,  d,  os  taken  tc  be  cr.e 
half  the  signal  wavelength  at  the  highest  frequency  (k-.in'-r*  "he  overlap¬ 
ping  signal  band,  with  total  length,  1  >>  The  accroach  cf  the 

study  is  to  isolate  interference  effects  resulting  fro.?,  the  closely 
spaced  configuration  of  sources,  and  to  do  this  in  a  setting  in  which 
reasonably  accurate  es-timates  cr  either  bearing  ray  be  expected  in  the 
absence  of  t.ce  interfering  source. 

In  order  to  assess  the  performance  less  resulting  from  the 
closely  spaced  source  configuration,  the  ratio  cf  Cramer-?, 
obtained  in  the  two  source  setting  to  the  bounds  o o ?:u t e d 
estimate  with  a  single  source  only  was  computed.  The  major  results  for 
the  situation  of  known  signal  levels  are  as  follows  : 

«r 

1.  Ever,  when  the  bearing  cf  one  of  the  sources  is  known  aprio- 
ri,  bearing  accuracy  for  the  source  at  unknown  bearing  is  strongly  degr¬ 
aded  if  the  SNR  of  the  known  source  exceed-  that  cf  the  unknown  source. 
The  degradation  factor  reduces  as  the  S'RR  of  the  known  source  increases 
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'.d  a  wall  de f i.ned  limit  and  then  saturates  at  a  nonzero 


both  source  bearings  are  unknown  and  the  two  signal 


noise  ratios  are  muon  s m a _ _ e r  than  t:.s  quantity 


then  _ the  accuracy  cf  either  bearing  estimate  is  degraded  by  the  fact; 
(Ij!"-  >>  1.  In  this  regime,  the  accuracy  of  either  source  bearing 

relatively  unaffected  by  the  SNR  cf  the  interfering  source,  provided  be 
signal  to  noise  ratios  are  less  than  t.ta  quantity  referred  tc  above*. 

3.  If  cr.e  of  the  signal  to  noise  ratios  exceeds 


\  1  " 


while  the  other  dees  net,  the  situation  is  quite  different.  In  this  case 
the  larger  SNR.  source  has  no  performance  degradations  resulting  from  the 
presence  of  the  second  nearby  source.  Hence  one  can  obtain  good  estima¬ 
tion  accuracy  even  in  the  presence  cf  interference  for  the  iarce  SNR 
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bearing  is  again  degraded  as  in 
signal  to  noise  ratio  "  regir.e,  with  cotimal  perfortancs  acheived  when 
the  power  ratio  is  one.  The  error  for  both  sources  is  proportional  to  the 
interference  source  power. 
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the  accuracy  reduction  over  the  known  signal  level  case  for  e : 
is  independent  of  signal  to  noise  ratio  and  equal  to  5.4  db . 
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6.  If  either  SNR.  is  larger  than  the  above  factor,  both  estim- 
degraded  over  the  single  source  accuracy,  unlike  the  result  det¬ 
ermined  earlier  where  the  larger  SNR  source  suffered  no  loss  in  bearing 
accuracy.  Incremental  errors  resulting  from  unknown  signal  levels  are 
proportional  to  the  ratio  of  the  larger  to  smaller  source  SNR. 
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the  lower  SNR  source  increases  with  the  interference  pow- 


7.  If  both  signal  to  noise  ratios  exceed  the  value  described 
earlier,  the  absolute  bearing  errors  rise  with  the “3/2  power  of  the  angu¬ 
lar  spacing  and  tend  to  zero  as  the  ratio  of  source  to  interference  power 
tends  to  infinity.  The  incremental  error  resulting  from  unknown  signal 
levels  increases  with  the  interference  tower  up  to  a  defined  value  and 
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:e  value 


9.  ?.e  achability 

sources  c ar.net  be 


cf  the  Cram.er-R.ac  lower  bounds  is  unlikely  if 
distinguished  (  i.  e.,  the  correct  number  cf 
sources  indeterminable  from  the  data  ).  This  question  was  pursued  by 
determining  the  probability  of  error  in  choosing  the  wrong  hypothesis  in 
the  following  binary  hypothesis  test  : 

Under  Hq ,  there  are  two  signals  present; 

Under  r>  ,  there  is  or.e  signal  present  radiating  the  total  racei- 
frem  bearing  angle  oq  located  at  the  weighted  average 
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Appendix  1 


Derivation  of  equations  (13)  -  (15) 


The  elements  of  the  Fisher  Information  Matrix  can  be  compute 
from,  the  following  expression  : 
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and  *  are  given  by  : 
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With  the  use  of  the  above  expressions,  one  can  show  that  : 
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Using  analogous  computations  for  element  J22  : 
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it  then  follows  that: 
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J22  is  identical  to  element  Jn  with  the  exchange  of  subscrip 
with  (2)  in  (A4)  : 


J22  ~  ^  Gk  |  P2S2  ( uk)  /N  Ufc)  |  ^2  “  Pl^2  f  +  pl2p22GkT2*i  i 


k  =  1 


Equations  (A3)  -  (A5)  are  equations  (13)  -  (15) 


Appendix  2 


-J 


Derivation  of  equations  (25)  -  (27) 


According  to  equation  (16): 

3  I  *  |2 

U  =  1 21*22 ! 


With 


N 

jvitV2|2=  ^  cos(u]{d(i  -  j)  (cos  (ax) -cos  (02)  ) /c) 


(A2.1) 


it  follows  that 


i ,  j  =  -N 


a~  iZl'tZ2l2=  ^  sin  (ax)  sin  (u^d  (i  -  j  )  ( cos  ( ax ) -cos  ( ax )  )  /c ) 

i,j  =  -N 

(A2.2) 

Introducing  the  notation  A  a  =  (ax  ~  <*2),  a  -  (ax  +  ax)/2,  along  with 
the  assumption  that  the  argument  of  (A2.1)  varies  negligibly  ever  the 
the  i  and  j  sums  one  obtains: 

!  2:1  "^2  i  2=  ^T((«k  d(i  -  j)/c}2-  (ukd(i3~  j)/c)4Aa2sin2(a)  +  .  .  )« 


With 


Aasin (a) sin(ax) 


.  .  2  4MN  (N  +  1)  (2N  +  1) 


(A2.2a) 


-  j  r  - 


i , j  =  -N 


For  large  M,  N  2  M/2  and 


2  .  M 
=  3  * 


i,  j  =  -N 


Truncating  the  series  of  (A2.2a)  after  a  single  term: 


1 21  t  22  I 


2  M2L2Aasin ( a) sin  ( ax )  2 


(A2.3) 


■  .*  «*_  ••  ••  •*.  -■  .■ 


The  expression  for  is  given  by  equation  (13): 
i’l  =  2Tr  \  (t-^—  (v*  v-  *)  )  (Z2V2i)  (  -r-^—  (vi  v  '  )  )  }• 

V  *  ^  J 

r  ± 

««  i  i  *  *  /  3  \  /  *  3 v i  .  »  *  ■  <  4  6  Vi  %  i  •  3  v  i  . 

=  2  i  ( v-  v«  )  ( v~< )  (v->  ■  - +  ( vi  v-> )  { vo  — — - )  ( v-?  -=— ) 

da 2,  ““  *•  da^  — «•  da^  da^ 

.  *  ,,Jvi  *  ,  ,Jvi  *  ,  ,  i  i2  3vi*  Jvi,  1  .  ,  _  . , 

(Z2  1  Zl )  Z2 )  t=^  Zl  +  V1  Z2  ^  (A2 . 4) 

d  a  ^  “  d  a 2.  “-*■  ■■*»  !  ^  d  ot^  f 

— £ , •  .  .  .  ~  d  V i  ijvi 

ror  far  fieio  sources,  v  vi  =  v  *>  v  i  .  A 1  so  — =•*■  v i  =  -  vi  — ==^-  . 

— *  — “  ~ *  —  dax  —  — L  dax 

Hence  the  middle  two  terms  of  equation  (A2.4)  are  identically  zero. 
Carrying  out  a  few  steps  of  algebra  the  remaining  terms  become  : 

N  -  o  2 


Ui  =  2  Y  u>: 
i  ,  j  ,  1  =  -N 


d  sin  (a-)  ,  .2  .  .  .  ,  .  1  (iasir.aduv) 

- ^2 (  1  +  1J  )  (  1  -  Jr - ^2 — ~ S— 


( i  -  j )  2  +  .  .  .  ) 


O  ^  n  -) 

M«—  <«  / 

at,  L  Sin  (  a  - 


+  0(/5.r2)2 


(A2.5) 


Finally 


I  .  =  -  — ^ —  -  |  V 1  ‘Vi  *  !  2 


d  a  x  <5  a  X 


ax  =  ax 


=  2  Kvi-viXfSl  *£1)  -  (21»»S1)  (*S1.  V.)[ 
^  ~  da  2,  da^  da^  da^  j 


2  2  2  N 

_  „  «v  d  Sin  V~  #  j  2  4J1  _  4MN (N  +  1)  (2M  +  1)  J2_._2,  , 

=  2  *~2 -  /  1  2  - 572 -  d  Sin  (^i) 


1, 2  =  ~N 


For  large  M  we  have 


,  _  M2L2sin2(ai)  2 

31  “  - cm2 - "k 


(A2.6) 


*Li _ I  v  1 1  vo  i  2  =  *k2d2sin(ax)sin(g2> 


j)2(  1 


(Aa Sin(a) (i  -  j) )  d  u>  ) 


)  +  ...) 


M2L2sin  (ai  )  sin  { a?)  , 

- ^  {  1  +  5  6  +  -  •  •  ) 


(A2.6) 


If  one  exchanges  the  index  1  with  the  index  1 
in  equations  (A2.3)  -  (A2.6),  expressions  for  Z2,  U2  and  *72  can  be 
obtained.  Substituting  the  variables 


xi  =  ( 


L2sin2{cn  )  .  1/2 


„  _  ,  L2sin2(a?) , 1/2 

x2  -  (  - 5c2  ) 


1  IZ. 'Z? 
f{2 


into  the  above  results  and  with  the  definition 
Kk  A  [  1  +  MS/N  +  H2S1S2/N2  I  ]_1 

substitution  yields  (25)  -  (27)  : 


rt  x 

Jll  =  ^  Hk  (MSi/N)  2X12«k2  [  1  +  0uk2Hk(MS2/N) 2 


(A2.5) 


k  =  1 


J12  =  M2SlS2/H2X1X2<1,k2  [  1  -  6 / 5 ^5 <wk 2  —  /5«k2HkM2S1S2/N2  ]  (A2.6) 


k  =  1 


X  n 

J22  =  Hk(MS2/N) 2X22«k2  [  1  +  ^«k2Hk(MS2/N)2  ] 


(A2.8) 


k  =  1 


Appendix  3 


vaticr.  of  equations  (39)  -  (41! 


Elements  of  the  Fisher  matrix  expressed  in  terms  of  integral 
12  a”d  I3  are  given  by  equations  (33)  -  (35)  : 


X12(MS1/N)2^  11  +  ^(MS2/N)2I 3  } 


( A3 . 1 ) 


r12  =  X1X2M2S1S2/N2  I  Ii  -  (6/5)  jJI2  -  /5M2S1S2/JI2I3  j 
r  2  2  =  X22(MS2/N)2  j  Ii  +  p(MS1/II)2l3  \ 


(A3 .2) 


A3 . 3) 


irect  integration 


II  = 


2  . 

M  CiJ 


1  +  MS/N  +  M*SiS2/N‘-  i  ) 


2  2 
/JW  /JW 

With  5  =  -  (  1  -  - 

2  5 


.  ) ,  terminating  the  Taylor  series  af te 


a  single  term  will  not  materially  affect  the  value  of  the  integral 
since  the  denominator  of  the  integrand  is  dominated  by  the  term  1 
MS/N.  An  approximation  for  I2  is  as  follows  : 


-  „  T  r  «  du 

1  ■  2  t  j  1  +  MS/N  +  M2S1S2/N: 


/!  W  ^  /  2  ) 


Changing  variables:  x  4  mk,  one  obtains  the  simpler  integral 


T_  1  -3/2 

2 1  1  +  MS/N 


I  n5* 


2  , 

x  dx 


XTT5I7H  “I ■—  f  1  -  (w^V^tarTW'?  } 


tar.  *  (W3X/  T  =  WX 


i  *  x  i 


i  w  x  ) 


One  then  obtains  for  : 


II  = 


15  =  -5- 


13  = 


2  7 

1  +  MS/N 

3  ( 

,r  ca 

Iculatior.s 

for  I2 

T 

1 

W2 

2  T 

1  +  MS/N 

3  K 

V21/ 

O 

*)  ) 

TW 

1 

w2 

2  7 

1  +  MS/N 

5 

T 

O  - 

1 

/  1  .i.  MC  /  V) 

_  w4 

"2  E 

3  2 

-  KWZ  +  .  .  .  ) 


2  - 1  2  - 1  / ' 

(  1  -  3  C.rx)  *(  1  -  (W“R)  ) 


{  1  - 


W  2r 


+  .  .  ) 


{  1  -  w2x  +  .  .  .  ) 


Substituting  (A3. 4}  -  (A3. 6}  into  (A3.1)  -  (A3. 3)  one  obtains 


jlt  =  TW  (MSx/N)2  W^  (  1  +  1_  kw2{  2S2  _ 

11  2a  1  +  MS/N  3  1  1  +  5  1  SX  1 


)  ) 


J12  = 

TW 

M  Si S2/N 

W 

<  1  -  is 

2  a 

1  +  MS/N 

3 

J22  = 

TW 

(MS2/N) 2 

W2 

(  1  +  ^  KW2( 

2  a 

1  +  MS/N 

3 

h*** 


2S- 

s2 


-  1  )  ) 


(A3. 4) 

'tar.  x  ( 

(A3. 5} 

(A3 . 6) 


which  are  equations  (39) 


(41)  . 


the  F.I.M  can  be  exp: 
defined  in  equations  (30) 


M  S-S?/N 


1  +  MS/N 


and  changing  variables  :  x  a  uR"' *  ,  write 


2  ,t  1  +  MS/M 


2  I  — 


x2dx _ 

-  0R~-x’/5  +  .  . 


(A4 .2) 


For  WR  <<  1,  terminating  the  5  series  after  a  single  term 
the  integral  value  only  marginally  since  the  denominator 
integrand  is  dominated  by  the  constant  term,  1.  When  W^l/2  > 
situation  is  potentially  different  because  the  integration  is 
over  a  region  in  which  powers  of  x  contribute  significantly 
denominator  of  the  integrand. 

The  integrand  of  (A4.2)  is  strictly  positive  and  mo 
since  the  integral  is  taken  over  the  positive  line  segment  ( 

)  the  integral  is  also  strictly  increasing.  Let  g(Aa)  ■  {  : 


ar  r  acts 
of  the 
>  1,  the 
carried 
to  the 

notonic ; 
0  , 


(M-l) /2 

g  {La)  =  1  -  |yj_*V2l^/M^  =  1  -  costujjd  -  j ) dias in ( v ) /c ) 

i.j  =  -  (M-l)/2 

A  Taylor  series  for  g  as  function  of  ia  is  given  by  : 


g{La)  =  g(0)  +  la  J  (0)  +■  La  g  (0)  / 2  L  +  .  .  . 


Because  the  La  -  dependent  part  of  g  is  a  sum  of  trigonometric  funct¬ 
ions  it  is  straightforward  to  demonstrate  that  the  following  proper¬ 
ties  of  the  infinite  series  hold:  The  odd  ordered  derivatives  of  g 


ft**.**  V mjCa 


M>  U) 


..  .  w  ..  3  _ 

U  .  'T 


e  r  i  v  a 
on ,  w 


=  C,  g  (C)  >  0, 

g  (0)  =  0 

r.ite  series 
■r.inating  th 

r.  alternating  inf i 
rrcr  caused  by  ter 

s.  The  error  can 

be  no  large 

the  (  n  +  1  )  terr 

:  of  the  ser 

as  a  result  bound 

the  i n  - 3 ^ i* 3 

s  strictly  positiv 

■e  and  nor.de 

bounded  as  well. 

Write 

g  ( 0 )  <  0 

an  ucoer  bound 


The  fact 


A 

a  ^  a**  * 

3  j  p  at 


;  if  {  =  ^-  +  £ 


sc  that  M“S1S2/hT^/ C  +  MS/M)  *  =  Xu"  +  Xc/{0/2) 
i  X  £  /  (  0  /  2  )  !  <  X4  X  ~  “  0  /  o  <  W4X  0 / 5 . 


With  X  =  u  X  “ '  " , 


Therefore  5  and  therefore  the  denominator  of  the  integrand  of  (A.  4 
can  be  upper  and  lower  bounded  : 


24  2-14 

1  +  x  +  pWX/5  >  1  +  x  +  X  fix  / 5  + 


2  4 

>  1  +  x  -/5WX/5 


which  yields  bounds  for  the  integral 


2  . 

_ x  dx _ 

1  +  x2  -  /JW4X/5 


2  , 

x  dx 

2  ~=^  4 

1+x  -  X  fix  /5 


2  , 

_ x  dx _ 

1  +  X2  +  0W4X/5 


It  will  be  necessary  to  consider  the  magnitude  of  /JW  x  /5  relative 


c:0)  (  2  t )  / 5  S3  tha- 


oduct  jJW  X/5  will  be  snails 
o  be  as  large  as  ICO,  the  signal 
be  as  laroe  as  3C  db  if  /5W-  is  to  be  sr.aller  th; 

iderations  we  will  assure  that  0WV\, 


=  WR*'*  -  (  1  -  W’jJ*  / 5 )  tan  ~  (W3- 


1  +  x*  -  W'R  V5 


The  lower  bound  is  given  by  : 


[  T  *  MS1/2  -  (  1  +  JSW4R/5  )1/2tan‘1(  Sr  K2 


Substituting  the  above  results  into  (A4.3)  one  obtains  for  the  upper 
bound  on  I]_  : 


_  r.ax  _  TW  1 

"  27  1  +  MS/N  { 


(W7?1/2)_1(l  -  W47?/J/5)1/2tan_1  (WK1/2)j 


I1ain  =  r  /MS/N  -  X_1  (  1  -  (W7?1/2) _1  (1  +  W4X/J/5)1/2tan~1(WK1/2)| 


x3  =  57  TT^ 


MS/N) 


X  (1  -  3a/4 (WX~' “) 


3«JJW3X1/? 

40  ' 

2 

The  upper  and  lower  bounds  for  1 2  and  1 3  differ  by  0 ( 0W  )  and  hence 
he  upper  and  lower  bounds  nearly  equivalent:  results. 


r  - 

12  "  2i  1  + 


.-1  1 


MS/N 


(A4.4) 


T  „  TW 
X3  =  07 


2,7  (1  +  MS/N)  - 


X~2  (  1  -  3.7/4  (W/?1/2)_1 


2  -1/2 

Let  y  4  ( W  7?  )  .Substituting  equations  (A4.3! 
(355  one  obtains  : 


-  {A4.55  into  (33) 


_  TW  2  (MSi /N) 


1  +  MS/N 


x"1  (  1  ♦  2|2- 


<  1  -  7fi>> 

(A4.6) 


*1*2 


m2SiS2/n2  -1,  . 


1  +  MS/N 


x  *(  1  -  y7  ) 


(A4.7) 


TW  „  2  (MS7./N)‘ 

2a  2  1  +  MS/N 


X-1  (  1  +  2|1 

s2 


y„/2U  +  §|i  )) 


(A4.8) 


Dropping  the  0(y)  terms  in  (A4.S5  -  (A4.8)  yields  equations  (48) 
(50)  . 


Appendix  5 


Derivation  of  equations  (55)  -  (71) 

In  this  appendix,  the  elements  of  matrices  J2  and  J3 
computed  using  the  following  expressions  for  arbitrary  elements 
the  Fisher  Information  Matrix  : 


VIT 

j.  .  —  |  y~ l  —  \  I  ”1  3Kv  ~r  “1  aKk  1 

Ji:  T‘  [  K  aI7  K  ae“  >  ~  )  T-  1  K<  K:<  “TeT  f 

k  =  1 

--  f-{ 


k  =  1 


From,  section  2,  Kk  and  K>  are  given  by  : 


Kk  =  Nk  I  +  Si_viv-  "  +  S2H2Z2''' 


_1  2  r  r  j. 

hk  =  jyr  T  ~  Gk  j  Pivivi  +  P2Z2Z2'  “  P1P2  (  ZlZl  *y’2v’? 


Z.2Z.2  *ZiZl  * )  f-  / 


where 


_  _  MS  1  /N  _  _  MS 2 /N  _  ,  _  _  ,  *  ,2  l“* 

?1  1  +  MSy/N  '  P2  1  +  MS2/N  '  °k  1  Pl?2  !-l  Z2  ! 


Calculation  o 


*  -  { K*-1  H* }  = 


We  have 


d  Kv  + 

aSi 


3  a  *  2 

also  let  y^-7  (  vjVjf  )  ft  Rj  and  - — 7  |z.l*Z2i  =  >lj-  Then  with 


wa  j  -  w  o  a  j 

=  Gk^  ?i?j  1j .  it  follows  that 


Jj—  =  ‘  PJ  {  Rj  ~  Pi  i— i f  +  ZiZi'Rj)  -  Pi7jGk(  Ziv^Pi 

-j-j*  pj  “  pipj ^ZiZi^ZjZj f  +  ZjZj fZiZi f ) ) j • 


Then  for  i  •  j 


=  — nj  (  1  -  2M?i  -  Gk?i(M  ?i  +  Pjjv^  v2  I  -  2?i?jMlv1  v2  I  ) 
Substituting  the  definitions  for  ?j_ ,  P  j  ,  £  and  Kk 

i  *  1 2 

t  -  Ivi  V-)| 

4  M- 


r*:<  (i  +  ms 2/:;)  (i  +  ms2/n)  .  1  +  MS/N  +  M  S; 

after  a  few  steps  of  algebra 

- Tr  1  Hr  Hf-1} =  ~2W—  &  ( 1  +  MSi/N  5 


For  i  =  j ,  =  vjVj4,  (A5.2)  becor.es 


1  +  MS/N  +  M2S1S2/N2^  j  s  Kk  1 


(A5.3 


|  dK'^  d  K>  l  _  Gk  P 
‘  l  as2  aaj  J  -  Nk 


«  “  /  O  i  u/T  —  |  *  _  —  .  *_  J* 

XTT  T  ■■•?-■  r  =  — Tr  i  Zi  vi  R i  -  Pi  (Ziv-  R-Z'Z' 
dSj  aaj  J  Nk  *w  -J-J  J  -  -« - - 

ZjZj  *ZiZi  *R  j  )  +  ^iijGk  ^°iZjZj  *ZiZi  *  +  NRjZjZj"”  ~ 

■> 

pipjZjZj‘ (ZiZ^ZjZj*  +  ZjZj"ziZir)  j 


(A5.4) 


Since 


Tr  |  ZjZj*Rj  j  =  0  and  Tr  j  ZjZj f (R jZiZi *  +  ZiZi^Rj)}  =  M  nj 


(A5.4)  becor.es 


nj  4  “MPi  +  P  jGk  (MP  i  j  v:  *v2  j  2  +  M2Pj  -  2M?iPj  |  vi_*v2  |^)  \ 


*  i  2 ,  I 


Nk 

so  that 


-«  {  $  f=j  <  ^‘l»  }  -  -  m-^n3 

Calculation  of  Tr  |  ||^  |I_(Kk"1)J  : 


(  1  +  MSi/N  5  )  n j. 


(A5.5) 


Calculation  oJ 


A  -ii‘^  -ni 


asj  k  "  Nk  as j  [ 


j  vjvj  -  Pj_{  viZi  v-iVj  +■  VjVj  ZiZi 


3Gk  1  i  _  * 

-ST"  P-'V^vj  +  P'V-v- 


-  ?^P j (ztZi  Z'Z- 


3?o  _  1 

3Sj  (1  +  MSj/N) 2 

=  -  G:</Ny  f 

(1  +  MSj/N)  -  !_ 


1 

Z-Zj^ZiZi")  [• . 

j 

a Gk  _  aG^  a?H  _  Gk  ivi_  vp !  p, 
aSj  -  a?j  TiTJ  “  (i  +  MSj/N) - 


[  A  5 . 5 


^  VjVj*(l  +  Pj?iGk  j  vi/v2  j  )  +  ZiZi  ’?i“  jZl 


I  2_ 

vo  |  Gv 


i  ■  ^ 

-?i  (ZiZitZjZjt  +  ZjZj  *ZiZi’)  d  +  Gkpi?  j  i  Zl  *Z2  !  )  [■ 


".'ich  1  +  P  j?iGk  |Zi"Z2  i  =  Gk'  ( A5 . 6 )  becor.es 


Gv2/Nv 


V'/  /*JV  j 


A  2  .  2  4.  A 

vizi  ^Pi  jZl  ”Z2  :  “  ?i  (ZiZi  ""ZjZj  X  + 


vjv,-  ~ZiZi  * )  S' 


c*  f=ll0“s  <*«  {  if 


■■  \  m  .  -<1^u2  ( i  ♦  MSi/K  e  i2 

L  dSj  dSj  ‘  J 


(A5.7) 


( A5  .  S 


Equations  (65)  -  (72)  can  be  obtained  from  equations  (A5.3),  (A5.5), 

(A5.7)  and  (A5.8). 


Appendix  6 


2 3  (  1  +  MS/N  ) 


After  direct  integration 


J  (  1  +  }<■ 

0 


I4  = 


2  t  (  1  +  MS/N 


(  -  T-TTP 


(w 2x  )  1/2« 


t  _  TW  1 

*5  2a  2(1+  MS/N  ) 


tan-1  (W  K1/2)  j 

I  _ L_ 

)*  \  1  +  W 


(w2*  )-1/2  * 


tan  1  (W  K1/2)  [ 


Fror.  Appendix  4 


_  TW 

"  2T 


1  +  MS/N 


r2  ( !  -  J-  (  w2r)'1/2x 


tar.-1(»s1/2)  ♦  2--x4w  } 


For  KW  <<  1  one  obtains  the  Taylor  expansions  : 


13  = 


5(1+  MS/N  ) 


(  1  -  W2R  +  . 


.  )  (A6.1) 


1 4  * 


3(1+  MS/N 


(  1  -  |  W2R  +  .  .  .  '  )  (A6.2) 


15  * 


(  1  +  MS/N  ) 


(  1  -  |  W2K  +  |  W4X2  -  .  .  ) 

J  D 

(A6.3) 


Substituting  (AS.l)  -  (AS. 3)  into  (65)  -  (71)  : 


M  SiS-j/N 

;i  +  ms7n) 1 

2M2Si S?/N2 
1  +  MS/M 


xi  r/2  -jj-  \  i 


MSo/N  - 


(A6 .4) 


TW  MSi/H  (  1  +  MSi/N  )  W2  „  .1/2  M  f  ,  6  ,„2 

27  - *(  1  7  MS/N  '  T2 - 3"  Xl  *  —  1  1  '  25“  M  " 


3/5W2  M2SiS2/N2  ", 

5  1  +  MS/N  ! 


(AS. 5) 


TW  MS  ->/N  (  1  +  MS2/N  )  W  v  1/2 

27 - "*  (  I  +'  ms7n~  F2 - 3~  *2  * 

3/5W2  M2SiS'>/N2  ", 


1  +  MS/N  j 


TW2  M2SiS2/N2 
2 t  (  1  +  MS/N  ) 


X2  /» 


1/2  M 


\  1  - 


/5W2  - 


(AS. 6) 


3/5W2 


MSX/N  - 


2M2Si So/N2 
1  +  MS/N 


6  /5W2 


TW  (  M/N  ) 

2 i  (  1  +  MS/N  ) 2 


f  ,  .  MS2/N  (  1  + 

\  3  1  +  MS/I 


MS2/N  ) 


(MS2/N)  2  (/5W2)  (1  +  MS2/N) (1  +  MS2/N  -  2MS]_/N) 

20  (1  +  MS/N) 2 


(AS. 7) 


TW  (  M/N  } 

2a  (  1  +  MS/N  ) 


{2  -  £ 


2M2SlS2/H2  ,  + 

1  +  MS/N  ' 


(^W2)2  ,  2  3(M2SiS2/N2)2 

20  1  5  (  1  +*MS/N  )2 


2M2SiS2/N2 


1  +  MS/N 


(  M/N  )4 
1  +  MS/N 


f  ,  JW2  MSt/N  (  1  + 

\  3  1  +  MS/N 


(AS .9! 


MSi/N  ) 


N)2  ^y2)2  (  U  +  MS!/N)  (  1  +  MSj/N  -  2MS2/N  ) _ )_)  1  (AS. 105 

m  on  '  /  1  +  MS/N  ' 2  '  • 


-  TW  (  M/N  ) 

-  l  2?  1  Tit  MS/fl 


t  t-  (  1  t  MS/M 


61  (flVQ 


O  O 

(fivt  “ 


f  (MS/N) 2  (1  +  MS/N)2*  2M2Si  S->/N2(l  -  2MS/N  -  3(MS/N)2] 
L  (1  +  MS /}!)*- 


Calculation  of  Xi^  : 

Xi  i  =  - -13 


(AS . 11 


2  -  _  _  o 

_a13  u  44 _ 2Ji3JijJ34  +  J14J37 

dec ( J3 ) 


2  -1 


consider  the  case  in  which  MS^/N  <  <  ( /jw  )  ;  i  =  1,2 

Atter  a  tew  steps  of  algebra,  one  can  shew  that 


:,A„  ,  f  W  1*  [  K_  1“  (  T  2  Vi4  r  .  ^  3  ,.2 

“3  L  27  N  (  1  +  MS/N  )  o  P  g  1  1  ■  5  £*  x 

MWv  -  2M2S1S2/W2)  1  _  12  flr.2 ^  pw2  (  MSi/N  )(  1  +  MS  1  /N  )  1 

2  1  +  MS/N  J  25  +  "3 - - - r~MS/N  ~ -  ^ 

J 

■Ti  i,’i  <  .T 1  «  =  -  T/Ll2  [  M  |  2  2  {  MSi/N  )2  MS2/N  (  1  *  MSi/N  ) 

2j  ,  l  N  J  (“I  +  MS/N  )b - - - -  * 


4  I  3  -> 

—  I  1  + 


{  1  t  ,3  ,„2  [  MS2/M  -  ]  .  ||  ,,2. 


-  2  m2Sis2/n2  I  1 

61  1  +  MS/N  J 


J  2j  =  r  TW  ]  r  M  (  MSi/N  )  (1  +  MS  1 /N  )  2  .  W 

14  33  L  2  ir  j  t  N  J  1  1  +  MS/N  )  6  X1  15  9~  * 


n  -  H  «W2  +  ^  MS2/N  (  1  +  MSp/N  )  _  6/5W2  M2SiSo/N- 

l  25  3  1  +  MS/N  5  1  +  MS/N 


,2  , 


After  considerable  algebra  one  can  show  that 


.  v  i-  %  , 


MS^/N  <<  [fiV  ) 


and  MS  2/N  >>  (AW  ) 


so  that  K/JW  <<  1- 


According  to  the  above,  MS 2 /N  >>  1  and  S2  >>  S]_  .  Based  on  the  res¬ 
ults  of  Appendix  5,  the  following  expressions  are  approximations  for 
elements  of  the  J2  and  J3  matrices  under  the  specified  conditions. 

2  3/2 

J13  »  -  ll  *1  -3-  B  KSi/M  (  1  -  2  fl  *  (MS 2/NJmV1  + 

+  |  (MS2/N)-1  ) 


:w  M*Si/N"(l  +  MSi/N)  1/2 
(MS2/N)*  * 


(  1  -  |t  A’.v2  -  |  MS;;./: 


I  ) 


T  ,  TV;  m  w“  -1/2  ,  . 

J23  »  27  N  1~  X2  P  (  1 


-  6  a';2  3  m 

aw  -  t  m 


25  w  -  J  MS i_/N AW  ) 


T  .  TW  M  S-  VT  1/2  ,  -  3  rr2  6  fffT2  , 

J2  4  3  2.t  N  So  3  *  (  1  10  lwl/J  ^  '  2:  ^  5 


2 

TV  M 


▼  -  a  n  kA 

=  27  N 


I|p-2(  1  ♦  §2  ,  ks2/h  ,w2  I’1  ) 


2 

TW  M  —2 

J34  *  U  (ms2/n) 


=  j44 


when 


It  follows  from  the  above  that  detJ3  =  J33J44-  with  J33  >>  J44 


MS2/N  >> 


(A*2) -1 


2  -1 
( A*  ) 


We  assume  that  the  above  statement  holds,  that  MS2/NAW  may  be  arbit¬ 
rarily  large.  From  section  3  : 


Xn  = 


2  2 
J13  J44  "  2j13j14j34  +  <*14  J33 


det  J3 


inserting  the  amove  approximations  vie.as  the  following  result. 

2  2 


=  ™  (MS:/M)2-^ - X-2{1  -  4  |i-  +  1-(MS2/:*0W2)“1-  f  (MSo/l 

a.  T  —  3  J  S'-'  J  3  “ 

^  ) 


-1 

N)  ‘ 


: r  aii  : 


v.  3~  44^2  3-  ~  ^14^34^23  ~  J3  4  Jt  4  w-_  3  +  Ji^J24j33 

det  J3 


~  ^  1  3  ^  2  3  _  £ 1 4 £ 2  3  _  J 2 4^1 3  +  4^24 

J33  J33  J33  J44 


tv  tj4*  c  T  -J  t  <r, 

=  ^  2MS-./M  X-X2  ~  (  1  -  4-  U'.r  -  i  MS./:i  /jvr  -  4 

2  s  -  -  *  3  2  o  5  -  St 


2  2 

Xtt  =  ?21  lAA  ~  1^.2.2J..2A'2i.  +...£Z1..  J21 

det J3 

2  2 
2  J 2 3  _  2  J2  3J2  4  +  ^24 

J33  J33  J44 


„  TW  20  ,-,.2.-1  Y  2  ,  12  _„2  6  /M  ,.,2 

s  27  y  (/J*  )  x2  — y-  '  1  "  t-=  Aw  ~  5  MSi_/N  0W  +  . 


Section  2,  equations  ( 3 S )  -{40)  give  the  elements  of  J2  . 
2-1  2-1 

For  MS2/N  >>  ( fiW  )  and  MSi_/N  <<  (fiW  )  these  become  : 


TW 

(MSi  /N)  2 

Xl2 

w2 

(  1  + 

2  T 

ms2/n 

*1 

3 

TW 

2  i 

MS 2/N  X]_X2 

2w 

3 

(  1 

6 

2  5 

TW 

2s 

MS2/N  X22 

W2 

3 

(  1 

♦§ 

ing 

\  terms  of 

inverse 

powe: 

t  MSt/N  /JW  -  t_  MS]_ /N  0W  ) 

o  q  ^ 

■  *  _  -7  MC.  /m  > 


,.,2 


TW  „  2  W  (MSi/N) 
"2  .  *1  *5  M  C  2  /  M 


fx  MSi/!I  0W“  ) 


J11  "  Xii 


{  1  + 


Xi  -> 


—  MS! /'I  (  1 


—  x2x2  - 

2  t  x  *• 


w*  +  ~  ms  2/u  yr 


25 

O 


tw  ^  y  i 

J22  -  x22  s  57  x2^  -3-  MS2/N  (  1  -  ^  /JW'-MS!/:;  ) 


At  lew  signal  to  noise  ratios  for  source  1,  i.e.,  MS2/N  <<  2 

variance  of  bearing  errors  are  given  by  : 


Var(ai )  5  (  57  X1 


75 


■N  ""  1 

jJW^  (MSj/N)  |i.  l 


S  2 

4*  J 


Var(S2)  >  (  57  X22  y|  W2  /5W2  MS2/N  } 


-1 


At  large  signal  to  noise  ratios  where  MSi/N  >>  1 


Var(51)  >  (  xx2  W2  ^(MSj/N)  ^  f 
Var  (s2)  >  j  Jr  x22  w2  ^5“  (  M2SiS2/N2)  j- 


2  Si 


i-l 


Equations  (A7.1)  -  (A7.4)  are  equations  (94),  (95),  (97) 


,  the 

(A7.1) 

(A7.2) 

(A7 . 3) 

(A7.4) 

,  (98) 


(1)  *u 


since  the  matrix  e 
terr.s  of  integrals  I3  -  i5 
approximations  for  the  inte 
A  method  for  obtaini 


egrals  . 


.  *e*  simpie  computations,  one  can  show  that  th° 

n.ust  lie  between  the  values  : 


i  +  ms/';,1  - 


k  *  i  wr 


2 t  (1  +  MS/M) 


1/2 

t/4  {  V?.  '  (  1 


!  V  ?  ' 

-  )  } 


per  and  lower  bounds  are  approximately  equal  because  we 
basis  of  physical  arguments  that 


2  1/2 
j5W  <<  W3 


Appendix  4  contains  details  of  the  argument. 

Cn  the  basis  of  these  approximations  and  the  expressior 
also  from  Appendix  4  one  obtains  the  following  ; 


It  = 


=  T*.  _ 1  -b~2  [  ,  _  3y t  ,  ,  .  /jw 

2  s  Tl  +  MS/N)  -  *  |  1  ~4 -  (  1  *  — 


—  _ L  .  R~2  (  i  _  3y*  j 

o  »  i  m  c  /mi  2  *  i  . —  r 


13  =  2  it  ■  (MS/N)  2  *  { 


*  1  1 
0-  1  ' 


Ji  3  = 


_  TW  -1/2  M  f  2 


2l  f  Xl  N  (  MS x /N  +  y  2  (  ~MS/N 


j,,  =  TW  -1/2  M  f  «  Si  1 

2l  1  N  {  Y  2  S2  MS/N  j 


3 

MS3/N 


=  -  - 1  /  2  _  /  *  s  ->  1  1 

2*  2  N  1  y  2  S L  MS/N  j 


-  -  •  _*•  *  id 

1 


\ 


w  <11 


tJjH  MH  c  m|H 


J  3  2  = 


:t  (MSi/n)*-  |  1 


4  (MS /:i) 


ms2/n  (  4MS2/N  +  2  MS2/N  )  r 


T  _  TW  (M/N)  yt 

"34  "  2t  (MS/N) -  4 


TW  (M/N) *  ;  , 

44  "  2 i  (MS2/N)*  t 


4  (MS/N)  2MS1/N  <  4MS2/N  +  3  MS2/N  )  }- 


:e:::3r.  of  the  above  approximations  for  the  elements  of  ^2  : 

2 

3 )  =  ^33^44  +  0 (y  ) 

2  2 
J  2  3 _  2J13J24J34  J  2  4 

J23  J33J44  J44 


7  * 


4  -  yi  (  3  + 


'14J34J23  _  J34J24J13  +  J14J2- 
J33J44  J33J44  J44 


(A3 .i; 


V  - 1 

T  t  *  *1*2 


y-t  (  1 


S*  So 

T2* 


-  /5  1  Xi2  *  4  -  y»  (  3  +  §-2- 

1  *  S  *■ 


3ns  (A3.1)  -  (A3. 3)  are  equations  (91)  -  (93). 


(AS .2) 


(AS . 3) 


1^2 

—  izi*v2'  '  one  obtains 

. •  —  ■ 


\  -  \  MS  +  M  S1S2  5i2  =  0 


The  following 

properties 

of 

the  solutions  to 

quadratic  equations 

yield 

2 

f  \ 2  =  MS, 

\l\2  = 

M 

sls2  *12 

M-2 

M-2 

2 

so  that  det (Kq ) 

11 

r' 

.  * 

+ 

K1 

*-> 

+ 

.to 

II 

+  MS  +  M  S]_S2  £12^ 

M 

2 

2 

=  N  Hk  ;  Hk  =  1  +  MS/N  +  M  Si.S2/N 

One  car.  use  ar.  identical  approach  to  compute  the  determinant  cf  the 
r.dtriA 


s5C- 


t .  \  ***  1  -j. 

(l-s)no  =  m7 

N; 


-  s  PqZOZO*  -  (1  -  (?iy2yj_ 


.  1 


P2Z2Z2  “  ?1?2 'ZlZl  ~Z2Z2  +  Z2Z2  "ZlZl  * ) 


S/M 

where  Pq  i  - . 

1  +  MS/N 

The  above  matrix  in  brackets  has  the  form  I  -  A.  The  matrix  A  has 
three  nonzero  eigenvalues  which  correspond  to  eigenvectors  of  the 
form 


avQ  +  bZl  +-Z2- 

The  matrix  I  -  A  has  (M-3)  eigenvalues  of  1  and  three  others  of  the 
form 


1  »  i  -  1,2,3 


where  the  are  the  nonzero  eigenvalues  of  A. 

Rearranging  the  terms  in  A: 

A  *  PoZOZO*  +  Gk?lZl(Zi*  “  p2 (Zl*Z2)Z2*)  +  GkZ2^Z2*  ~  Pl(Z2*Zl)Zl  ) 

and  writing  the  matrix  equation  det(A-M)  =  0  one  obtains  the 

following  expression.  For  brevity  of  notation,  define 

Pi j  A  Zi*Zj  for  i >  3  =  0,1,2. 


tr  w 


cons.ueraDie  a^genra  one  can  show  that  the  above 
equation,  results  in  the  following  cubic  equation  in  \  : 

*3-  K.  (i  +1MS/:n  [  (1  +  MS/N) {MS/N  +  2m2s1s2/n2512)  -  sm2s1s2/n2* 


(2  4-  MS/m 


K^nV^s/N)  (  11  '  s)  11  +  ::s/n)m2siS ;/:>2512  - 


sMS/N  (  MS  i_/N  501  +  MS2/N  502  +  M2S1S2/:J2(2  §12  -4  501  +  £C2 


h,(i!;  ms/'i)  “s/"  k2*;*:/»2  [  512  +  SOI  *  t< 


-  2  X 


U9.1) 


'allowing  notation  is  used  in  (A9. 


~  |Piji  ar*d  X  =  1  -  —  Re  (P01P12P2Q^ 


:e  only  the  suns  and  products  of  the  Xj_  enter  into  the  determinant 
not  necessary  to  solve  for  the  individual  eigenvalues.  Write 
ove  equation  as  follows  : 

3  2 

X  +  aQ  x  +  a2  X  +  a2  =  0. 


X l  +  x 2  +  x  3  =  -  ag;  X  i  X  2  +  x  2  X  3  +  X  2  X  3  =  a  2  ;  X i X 2 X 3  =  ~  a  2 


(1  +  MS /II)  Hk  [ 


1  +  MS/I!  ->■  s  M  S2.S2/.I  Si; 


;)  MS/N  (  MSi/II  S01  +  ms2/n  S02  )  +  s  C-s)  (  $12  +  50i 


h  the  previous  results  for  det(Ki)  and  det{K0),  the 
;i  ( s  )  fcecor.es  : 


in  [  i  +  K  S 12  J  " 


3^12  -  s(l-s)  (  Sol  - 


)  R  +  s“(l-s)  MS/N  R  U12  +  S02  +  Soi  “  2  X)  (A9.2) 


is  equation  (107; 


New  consider  equation  (A9.2)  with  the  specialization  of  vq  : 


~  0  a  ;<  v  1  _  0  u>  ;<  v  2 

ZO  =  C  1*  e  ,  e 


“0  6J:<VM-1  T 
-  e  ] 


sial  +  s2a2 

=  d  s in  (  -  ) /c 

S 


the  results  of  section  2  : 


y 


cos  (  ^vd  (  i- j  )  ( cos  an  -  cos <32  )  /c ) 


addition,  WLsir.(a)  lct/c  <<  2t,  one  obtains 


,  _  V  1  ,  2,2  .  2,  ,  2,  ,  .  .,2 

512  ~  /  2  d  sm  ^  Aa  ^ 

i  j  =  Jill 

i,  J,  2 


2  2  2  2  2 
=  L  sin  (a)  4a  /5c 


lOl  =  1  -  1^-  ^  cos  {  (cosao  -  cosen  )  { i-j  ) /c  5 

•M.  M-l 
1,2  =  -2— 

Sl~S2 

Using  the  previous  definitions,  ao  =  a  +  A  -  and  it  therefore 

S 

follows  that 

si-s2 

cosaQ  -  cosai_  =  cos(a  +  h  -  )  -  costa  +  h)  = 

S 

S1-S2  s2  s2 

=  h ( 1  -  - )sina  =  2h  —  sina  =  4a  — 

s  s  s 


After  a  few  steps  of  algebra,  it  follows  that 

2  2  2  2  S2  2  2 

$01  3  6)^  L  sin  a  ha  {  —  }  /6c 

S 

2  2  2  2  Si  2  2 

Similarly  $02  s  ®k  L  sin  «  ( — )  /6c 

S 

2  2  2  2 

Let  pi  L  sin  a  ha  /6c  .  Then 


, '  »  "  ,  -  -  "  -  '"o  ■  *  •  ,■  ‘  '  •  "  _*■  *,*■  *  -  **•  ’•*  *  *  -”  «*  *•*/,*  V 

fiV  «  J  O' -  ■  •  ■  -  "  *-  *  -  -  *  v*  *.  ■s,  S  l.%  '*  *«  .V  wVlA.  -N  -V  i  %  . 


»*riWfirr  it  rr  ^;iw:iVrrr»,r»,  < 


*12  3  P 


c'  nally  for  X  : 


*01  =  «k  0  Sy  /s 


X  -  1  -  —  Re  (  ?oi?l2  ^20  )  = 


o  no 


*02  =  0Jk  S1  /S 


AS  .3 


=  1  "  =o 


l - 


d  ' 

exp  (  -jufc  —  |  itcoSaQ  -  cosay)  +  j(C0Say  -  COS  ay)  + 


i  /  j  ,  1  = 


1 (cosay  ~  COSoq) |  ) 


c- 

=  1  "  S3  y  003  (  \  **  sin*  j  i  |2  +  j  "  1  f1  M 

4  ,  -  M-l 
-  2 

Under  conditions  in  which  LW  ia  sir.a  /c  <<  2s,  one  car.  expand 
cosine  function  in  a  power  series  resulting  in  : 


J  Uk2d2da2sir.2a  /  M3c2)  y  (i(S2/S)  +  j  -KSy/S)) 


i,j,l  = 


Because  the  i,j,l  suns  are  carried  over  the  synnetric  interval,  only 
terns  with  even  powers  contribute  the  suns.  Using  the  definition  for 
P  it  can  be  shown  that 


2  2 

1  2  S  y  S  2 

X  S  —  uk  p  (  1  +  —  +  —  ) 

2  S  S 


( A-9. 4 ) 


Returning  to  equation  (A9.2),  by  inserting  the  results  of  (A9.3)  and 
(A9.4)  one  obtains  the  sinplif ication  : 

*12  +  *01  +  *02  -  2X  =  0 

S  S  2 

—  *01  +  —  *02  = 

S2  Sy 


[•] 


Ms)  = 


j  In (  1  + 


1  +  X  fi  01  >  “  S 


For  large  TV  products,  one  can  approximate  the  k-sums  by  integrals. 
Ms)  =  ^  f  [lr.(  1  +  Zfi*2)  -  - ~  ■; -  diJ 

*  *  L  i  ^  /(puJ^S 


Changing  variables  : 


V  { 3  j5 ) 


W  (  3  0  S  } 


T  f  -1/2  f  2  3  -1/2  T  x  dx 

Ms)  =  — ',{?.?)  I  ln(l  +  x  )  dx  -  (Xfis  )  !  - 

2M  i  ]  1  +  x^ 


tv  r 


2  -1/2  -If  _  2  1/2 


=  —  i  lr.(l  +  XfiW  )  -  2(1  -  (XfiW  )  tar.  I  (KfiW  ) 


-1  2  1/2  -If  2  1/2  1  ] 

-  s  (  1  -  (XfiW  s)  tan  I  (XfiW  s)  )  j 

i  J  ) 

2  -1 
For  R0W  <<  1.  power  series  expansions  of  the  in { * )  and  tan  (*) 
yield  the  following  polynomial  in  s  : 


TV  r  (XfiW 

Ms)  =  —  !  -  - 

2 *  [  10 


2  2  2  3  2  3 

(XfiW  )  (XfiW  )  (R0W  )  2 

-  s  +  -  -  -  s  + 

5  21  7 


One  can  compute  the  root,  sr,  for  any  order  of  the  above 
polynomial.  However,  it  can  be  shown  by  straighf orward  calculations 
that  the  nfc"  term  in  the  s  polynomial  is  proportional  to 
2  n  +  2 
(XfiW  ) 


so  that  corrections  to  the  location  of  the  root  by  including  higher 
and  higher  order  terms  will  tend  to  vanish.  The  first  order 


